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F.  H.  Fenlon 


SUMMARY 


Following  a  review  of  the  approximate  models  previously  used 
to  provide  parametric  scaling  laws,  an  exact  solution  for  the  asymptotic 
far-field  pressure  of  a  parametric  array  is  derived  from  the  plane  wave 
form  of  Burgers'  equation.  The  influence  of  spherical  spreading  losses 
is  then  included  by  matching  the  solution  at  low  primary  wave  amplitudes 
with  a  spherical  wave  solution  obtained  by  the  method  of  successive 
approximation.  From  the  matched  asymptotic  solution,  the  referred 
pressure  of  the  difference-frequency  signal  at  one  meter  from  the  source 
(i.e.,  the  equivalent  difference-frequency  source  level)  is  derived  aB  a 
universal  function  of  scaled  primary  wave  amplitudes,  frequencies,  source 
dimensions,  and  physical  constants  of  the  medium.  With  the  aid,  this 
function,  which  defines  the  saturation  limit  for  a  parametric  array,  the 
maximum  conversion  efficiency  r»r.  be  evaluated  for  a  particular  set  of 
input  parameters.  The  amplitude  dependence  of  the  directivity  function 
and  directivity  index  for  "spreading  loss  limited"  parametric  arrays 
can  also  be  determined.  Scaled  performance  characteristics  obtained  from 
the  model  are  shown  to  be  in  good  agreement  with  data  published  in  the 
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INTRODUCTION 


order  to  facilitate  the  design  of  parametric  arrays,  different 
investigators'*'  ”*  guided  by  experimental  observation  have  constructed  per¬ 
formance  characteristics  from  simple  approximate  solutions  of  the  second- 
order  nonlinear  wave  equation.6  These  solutions  comprise  (1)  Westervelt ' s7 
far-field  approximation  for  the  difference-frequency  signal  generated  by 
nonlinear  interaction  of  infinitely  plane  uisaturated  primary  waves  of 
finite  amplitude  (i.e.,  waves  whose  peak  amplitudes  are  below  their 
respective  shock  thresholds),  subject  only  to  viscous  absorption,  and 
(2)  the  corresponding  approximation  for  unsaturated  spherical  primary 
waves  derived  by  Cary8  and  the  author.8,9  As  shown  in  Fig.  1,  the  primary 
wave  fields  of  a  finite-amplitude  source  of  area  Sq,  embedded  in  an  infinite 
rigid  baffle,  operating  simultaneously  at  angular  frequencies  and  can  be 
treated  as  plane  collimated  waves  within  a  distance  from  the  source 
rQ  =  sQ/*r  (de"  bribed  as  the  "collimation  distance"  or  "Rayleigh  distance") 
where  is  the  wavelength  of  the  mean  primary  wave  frequency.  At  distances 
greater  than  r  the  primary  fields  can  be  represented  as  spherically 
spreading  waves.  If  the  near  field  primary  wave  absorption  loss  o.  r  is 
large  enough  to  ensure  that  the  primary  waves  are  sufficiently  absorbed 
within  rQ  to  the  extent  that  no  further  nonlinear  interaction  occurs 
beyond  this  range,  the  parametric  array  is  said  to  be  "absorption  limited," 
and  Westervelt  s  solution  can  be  used  in  this  instance  to  determine  the 
pre-shock  asymptotic  far— field  form  of  the  difference— frequency  signal. 
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On  the  other  hand,  if  ot^r^  is  very  small,  the  primary  wave  interaction 

takes  place  primarily  beyond  r^  and  is  limited  essentially  by  spreading 

8  9 

losses,  so  that  in  this  instance  the  spherical  wave  solutions  ’  can 
be  used  to  define  the  pre-shock  difference  frequencv  pressure  in  the 
far-field  of  the  parametric  array.  Within  rQ  a  virtual-end-fire  array 
is  formed  whose  half-power  beamwidth  0  (r)  at  pre-shock  primary  wave 
amplitudes  decreases  with  range  as  Jr.  If  the  array  is  "absorption 
limited"  its  final  length  is  entirely  determined  by  the  primary  wave 
absorption  coefficient  a^,,  so  that  the  limiting  far-field  beamwidth 


0_  =  A_aT  =4/aT/?.k_. 

Alternatively,  if  nonlinear  interaction  also  takes  place 

beyond  rQ,  the  end-fire  arrav  length  is  curtailed  both  bv  viscous 

absorption  and  by  spreading  losses.  In  this  instance,  as  r  increases, 

the  half-power  beamwidth  0_(r)  asymptotically  approaches  that  of  the 

primary  wave  product  beam  pattern  0  formed  beyond  r  ,  where  6  m  u/k  a 

o  o  o  o 

for  a  circular  piston  source  of  radius  a  with  k^a  >  1,  operating  at 

pre-shock  primary  wave  amplitudes.  Eventually,  0_(r)  equals  0q  at  a 

range  r*  =  r  (w  /w  ) . 

0  o  o  - 

3 

Likewise  as  shown  bv  the  author  it  follows  that  if  a  r'  >>  1. 

T  o  ’ 

the  parametric  array  is  absorption  limited  so  that  the  far-field  half¬ 


power  beamwidth  is  9_  at  pre-shock  primary  wave  amplitudes.  Alternatively, 

under  the  same  conditions  if  a„r'  <<  1  then  0  will  be  the  far-field 

1  o  o 

beamwidth.  Thus,  by  identifying  as  the  key  parameter  which  determines 

whether  a  parame.  ic  array  is  "absorption  limited"  or  "spreading  loss 
limited",  Mellen  and  Moffett  ’  combined  the  simple  plane  and  spherical 
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wave  solutions  to  obtain  an  approximate  asymptotic  far--field  solution 

for  all  values  of  a_r'.  This  solution  can  be  expressed  as, 

1  o  1  ’ 

2 


(a  / 2)  {  /ui  )' 

1,2  o 


r'  .  -nr'  ,  , 

o  2,  T  dr'  .  ,  ,  . 

rN(r  )e  r - 

o  o 


„  -a  r  ,  r  -a  r 

Tp(r 1  )e  T  £-}<-£)«  -  , 

i  r  r  r 

r  (i) 


where  all  variables  are  defined  in  the  list  of  symbols  and  p  is  normalized 

with  respect  to  the  peak  amplitude  at  the  source  p^,  rather  than  to  either 

i  o 

of  the  individual  primary  wave  amplitudes,  as  in  the  Mellen  and  Moffett  ’ 
model  (where  the  latter  are  assumed  to  be  equal). 

If  the  amplitude  taper  functions  and  T  ,  which  account 
symbolically  for  finite-amplitude  absorption  of  the  primary  waves,  are  set 
equal  to  unity,  the  terms  in  Eq.  1  can  be  combined  and  integrated  in 
various  ways  giving  the  unsaturated  solutions  shown  in  Table  la.  Inspection 
of  Fig.  2  also  shows  that  these  unsaturated  solutions  are  in  good  agreement 
over  a  wide  range  of  o^.  Mellen  and  Moffett's  model1’2  goes  beyond  the 
unsaturated  case  however,  by  providing  an  exnlicit  expression  for  the  ampli¬ 
tude  taper  functions  and  Tp,  which  are  assumed  to  be  of  equal  weight. 

This  expression,  which  is  of  the  form,  TN  =  Tp  =  l/l+(o/2)2,  was  obtained 
bv  empirical  observation  of  the  extent  to  which  the  peak  amplitude  of  a 
monofrequency  finite-amplitude  wave  is  reduced  due  to  phase  advancement; 


the  "distortion  parameter"  a  being  determined  by  .e  extend  to  which  the 
peak  of  the  waveform  advances  in  an  inviscid  fluid, 

ft  d(r'/r_) 


rr 


[ i. e. ,  a  -  a 


o  dr ' 

J°o 

o  o 


rr 


dr'  „ 
ZT~  ~  O. 


=  a  Sinh_1(r/r  )]. 


/l+Cr'/r^2 


Substitution  of  such  a  taper  function  in  Eq.  1  however,  assun.es 
that  the  primary  fields  can  be  treated  as  independent  monofrequency  waves 
losing  energy  onlv  through  transfer  to  self-generated  second  harmonics. 


- - . .  ...  -  ■  - - 


,  — - —— —  • 


Ihus,  in  making  this  assumption  the  additional  loss  of  energy  transferred 
to  higher  self-generated  harmonics  is  neglected,  in  addition  to  degenera¬ 
tive  coupling  of  the  tvpe  discussed  by  Tjotta,10  ilobaek  and  Vistrheim,11 
due  to  interaction  of  the  sum  frequency  component  and  the  primary  waves 
themselves.  That  such  effects  are  not  inconsequential  is  shown  bv 

Q  I  o 

numerical  solutions  of  Burger's  equation,  and  bv  Blackstock's  use  of 
the  dual  frequencv  Bessel-Kuhini  series  (which  holds  in  an  inviscid 
or  weakly  viscous  fluid  prior  to  the  critical  range  at  which  shock  forma¬ 
tion  occurs)  to  demonstrate  the  "absorption  of  sound  by  sound,"  which 
is  a  direct  consequence  of  cross  coupling  between  the  primarv  wave  fields. 

In  a  subsequent  paper,  Merklinger,  Mellon,  and  Moffett1^*  derived 
a  more  sophisticated  taper  function  for  spherically  spreading  monofrequencv 
waves,  which  unlike  the  previous  function,  includes  the  effect  of  viscous 
absorption  on  the  phase  advance  of  the  peak  amplitude.  This  function  was 
obtained  from  an  ordinary  nonlinear  differential  equation  of  the  Bernoulli 
form,  wnich  Merklinger  '  originally  deduced  bv  ciuasi  analytical 
arguments  in  order  to  provide  a  finite-amplitude  taper  for  plane  mono¬ 
frequency  waves.  Actually,  Merklinger's  equation16, 17  had  previously 
been  deduced  by  Nestervelt18  and  solved  by  Wiener19  for  the  case  of  finite- 
amplitude  standing  waves.  Furthermore,  it  is  easily  shown  that  this 
equation  is  simply  a  reduced  form  of  Burger's  equation  in  the 
spectral  domain  for  a  monofrequency  wave,  when  all  self-generated  harmonics 
higher  than  the  second  are  neglected. 

Nevertheless,  substituting  the  taper  functions  shown  in  Table  lb 
(which  also  gives  their  integrated  forms)  in  Eq.  1  leads  to  the  approximate 
expressions  outlined  in  Table  Ic  for  the  asymptotic  far-field  pressure  of 
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a  "partially  saturated"  parametric  array;  the  so  called  "Modified  Berktay- 

Leahy  Model,"  derived  by  the  author,  being  included  for  the  sake  of 

completeness.  For  <<  1  these  expression  reduce,  as  required  to  those 

of  Table  la.  It  should  be  noted  in  this  connection  that  although  Bartram's 

model5  has  been  included  in  Table  Ic  iL  is  only  valid  for  "absorption- 

limited  arrays"  (a.  >>  1)  since  it  is  not  asymptotically  matched  for 

all  values  of  m_r  . 

T  o 

Referring  to  the  expressions  of  Table  lc,  the  term  "saturated" 

1  2 

was  previously  used  ’  to  denote  the  fact  that  p  (r)/pQ  becomes  constant 
for  oqA>>  1.  This  use  of  the  term  is  misleading  however,  because  even  if 
Pl(r)/P0  becomes  independent  of  pQ,  p_(r)  remains  linearly  dependent  on  p 
for  all  subsequent  values  of  a  A  (a  parameter  directly  proportional  to 
PQ)  so  tllat  the  svstem  is  in  no  way  "saturated."  Hence  the  term  "partially 
saturated"  has  been  used  instead  to  denote  the  linear  dependence  of 
p_(r)  on  p  .  Clearly,  the  continued  dependence  of  p'(r)  on  p  is  an 

w  -  O 

unsatisfactory  result,  since  it  violates  the  nature  of  stable  nondispersive 

nonlinear  systems  which  eventually  tend  to  approach  well  defined  states  of 

equilibrium,  as  demonstrated  for  example,  by  the  results  of  Shooter,  Muir, 

20 

and  Blackstock  for  the  case  of  monofrequency  finite-amplitude  waves. 

Thus  p^r)  must  eventually  approach  "saturation"  as  Pq  increases  indefinitely. 

2 1 

In  the  final  analysis,  as  Cary  has  pointed  out,  the  real 
weakness  of  the  solutions  outlined  in  Table  Ic  is  that  there  is  no 
criterion  for  determining  the  conditions  under  which  they  fail,  nor  is 
there  any  means  of  establishing  the  errors  involved  in  neglecting  all 
finite-amplitude  loss  mechanisms  other  than  second  harmonic  generation. 

10 
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I 

I 

[  For  these  reasons  another  approach  is  required.  Tl  lis  will  be  the 

subject  of  the  following  section  which  makes  use  of  Burgers'  equation 
1  to  determine  the  amplitude  dependence  of  the  asymptotic  far-field 

y  pressure  of  a  parametric  arrav,  and  relies  on  the  more  general  second- 

order  nonlinear  wave  equation  to  determine  its  frequency  response. 

Unlike  the  models  outlined  in  Table  Ic  however,  the  effect  of  all  possible 
spectral  interactions  on  the  far-field  solution  are  included.  Thus,  the 
major  assumptions  on  which  the  model  is  based  can  be  stated  as  follows: 

—  (1)  As  in  the  case  of  all  previous  models^  5  the  amplitude  and 

frequency  responses  of  the  asymptotic  far-field  solution  are  assumed  to 
be  independent.  This  implies  that  the  far-field  difference-frequency 

»  » 

pressure  p^(r)  can  be  expressed  symbolically  in  the  form, 

n  ~nL-r 

Pl(r)/Pa  =  F(o1  2,  aTro)  (rQ/r)e  ,  1  <  n  <  2. 

(2)  As  in  tlie  case  of  all  previous  models^  ^  a  nondispersive  medium 
is  also  assumed. 


1.  THEORETICAL  ANALYSIS 


As  shown  by  Blackstock  the  propagation  of  one  dimensional 

progressive  finite- amplitude  wave  in  nond ispersive  thermo  viscous 

fluids  can  be  described,  correct  to  second-order  terms,  by  means  of 

22 

Burgers'  equation.  Using  a  slightly  modified  form  of  Blackstock's 
stretched  coordinate  system,  this  equation  can  be  expressed  in  terms 
of  the  excess  pressure  p',  normalized  with  respect  to  its  peak  value 
at  the  source  pQ,  for  plane  cylindrical,  or  spherically  spreading  waves 
as , 

2 

pip  r)  P  V 

^  -  p  aTT  -  (l/A)  ~^2  -  0;  P  =  (r/ro)\P'/Po),  (A  -  0,  1/2,  1),  (2) 

3 1 

where  the  variables  C,  t'  (each  with  the  dimension  of  time)  and  the 
parameter  A  are  functionally  related  to  measurable  parameters  in 
Table  II;  the  value  of  9,  being  determined  by  the  type  of  wave  under 
consideration,  i.e.,  plane,  cylindrical,  or  spherical,  respectively. 

TABLE  II  -  Functional  dependence  of  £,  t'  and  A  on 

r/rQ  and  t,  for  plane,  cylindrical,  or 

spherical  waves;  o'  =  (Bp  r  /p  c3); 

A  =  a’/fir  .  °  0000 

o  o  o 

Wave  Type  9,  c/a  t'  A 

o 

Plane  0  r/r  t  -  (r/c  )  A 

o  o  o 

1/2  “1 
Cylindrical  1/2  2[(r/r  )  7  -1]  t  -  (r-r  )/c  A  (l+C/2a') 

o  o  o  o  o 

Spherical  1  In (r/r  )  t  -  (r-r  )/c  A  exp(-C/o') 

o  o  o  o  o 


F* 


■  \ 


I 


1 


I 


I 


I  | 
1 


a- 

•**. 


I 

»  » 


!  i 


1 


T 


I 


r 


Now  the  plane  wave  form  of  Eq.  2  can  be  cons iderably  simplified  by  means 

r  23  24 

of  the  Hopf-Cole  transl 

dependent  variable  ip  as, 


23  24 

of  the  Hopf-Cole  transformation,  which  expresses  P  in  terms  of  a  new 


P(C,f)  =  (2/A  )  ■—  {ln^(A,t') 

O  o  L 


ft1 


nr 


=  exp{  (A  /2) 


P(C,t")dt"l . 


(3a) 

(3b) 


Substituting  Eq.  3a  in  Eq.  2  results  in  the  linear  heat  conduction  equation, 

(4) 


-  (i/ a  )  ^L.  =  o. 

dC  o  ^,2 


at 1 


Taking  the  Fourier  Transform  of  Eq.  4,  and  solving  the  ordinary  differential 


equation  thus  formed  gives. 


2  / 

'v  %  r,/A 

*u»(°  =  V°)e 


where 


- 


.tOe'^'df 


and 


-  IT 


*  U)eJ“E'du. 
0) 


(5) 

(6a) 

(6b) 


If  Eq.  5  is  substituted  in  Eq.  6b,  Mc.t')  becomes, 

1 


^ (C, t')  = 


2tt 


% 


-u‘c/ A„ 

ip  (o )  e 

0) 


dm. 


(7) 


where  ip  (o)  is  given  by  Eas.  6a  and  3b  as. 


i(o)  = 


ip(o,t')e  dt 1 


“  (A  /2) 
e 


P(o,t")dt" 
o  e'ja,t 


dt'  . 


(8) 
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In  the  case  of  a  parametric  source  operating  simultaneously 


at  angular  frequencies  to^  and  the  normalized  excess  pressure  at  the 
source  is  given  by  the  expression, 


P(o,t)  =  P  sin(u  t  +  4»_ )  +  P  sin(to„t  +  $  ) ;  P  =  (p  /p  )(i  =  1,2).  (9) 


Substituting  Eq.  9  in  Eq.  8  and  expressing  the  resulting  exoonent  as  a 

25 


double  Bessel  function  series  of  imaginary  argument  gives, 

f  *  j  (to  -to)t'  .1  d>  {(1-/2)003  -J),  +  (ro/2)cos  v  } 
nm  ran  L  12  2 


*„<»)•  n 

n,m=-°° 


e 
—  00 


e  e 


2tt  y  T  (-l)n+ml  (r./2)i  ( r  /2)5(ai  -to)e~  '  ™e 
L  n  i  ra  2  nm 


^nm  ^r^/2)cos  4>  1  +  (l2/2)cos 


,  (10) 


n,m=-°° 


where  the  Dirac  delta  function  6  (to)  is  defined  as, 
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6  (to)  =  2^- 


<ii) 


with  tonm  =  (nto^  +  mto2);  $  =  (n^  +  m^);  n,m  =  0,  ±1,  ±2, 


(12) 


and 


ri  •  (V“i)poi  ‘  a -i.2).  <i3) 


If  Eq.  10  is  substituted  in  Eq.  7,  iK£,t*)  becomes, 


00  -a)  £/A  j  (u>  t*  +  6  ) 

<Ks,t')  =  l  l  (-1)  (r  /2)I  (r  /2)e  11111  0  e  nm  nm 

_  n  1  m  z 

n,m— ” 


{(r^Jcos  4>1  +  (f  2/2)cos  <{>2 } 


.  (14) 


Returning  to  Eq.  3a  and  using  Table  I  (with  i  =  o)  to  replace  t;, 
P(r,t')  becomes, 


^  4  anm^Un, 

P(r.t')  =  -(4/A  - 


-a  r  -a  r 

e  n*  '  sin  (to  t'  +6  )  +  to  e  n’  m  sin  (to  t '  +  4  )} 

+m  n,+m  n,+m  n,-m  n,-m  n,-m 


-a  r 


"  -a  r 

1+2  ll  a  {e  n’  cos  (to  t'  +  <j>  )  +  e  n»  ra  Cos(to  t'  +  d>  )} 

L  nm  n,+m  Tn,+m'  n,-m  vn,-m 

n,m=o  ’  ' 


(15) 


14 
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i 


1  i  |i 

I 


-  . - 


- 


n.  I  (r./2)T  (!'  /  2) 

,  ,  , . n+m  r  n  1  m2,  .  1 , . 

where  a  =  (-1)  1-= — 77 — -pp-z — tt— r^vd  •  U«) 

nm  T.  (I  / 2 )  I  (1  /2) 

o  1  o  2 

Inspection  of  Eq.  15  shows  that  each  term  in  the  double 
summations  of  the  numerator  and  denominator  decays  exponentially  at  a 
rate  dependent  on  its  absorption  coefficient.  Since  the  latter  is 


proportional  to  the  frequency  squared,  it  is  clear  that  as  r  increases 
indefinitely  the  term  with  the  smallest  absorption  coefficient  (which 


in  this  instance  is  that  of  the  difference  frequency  component  w.  . 

1  ,  1 

will  outlive  the  others).  Thus,  the  asymptotic  form  of  Eq.  15  becomes. 


P (r , t ' )  -  -  (4/Aq){ 


aq  sin (w  t'  +  <(>_) 


— } ,  a  r  >  1 


(17a) 


1  +  2a  id_e  cos(lo  t'  +  <J)_) 

1,1  ~ 


-2a  r 


-  ~  (4/AQ)a1  _1<Jj_{e  sin(u)_t'  +  <j>_)  -  -j  e  sin  2(w_t'  +  <{>_)  +  . 

(17b) 

where  and  a_  =  .  (18) 

At  still  greater  ranges,  only  the  first  term  of  the  distorted 
sine  wave  described  by  Eq.  17b  survives,  so  that  eventually. 


P(r,t')  -*  P  (r)sin(uj  t  +  tp  ) ,  a  r  >>  1 , 


where 


P  (r)  =  -(d)  /id  )  (4/r  )e  ;  F  =  (A  /a>  ) ;  u>  =  -r-(id..  +  io0) 

-  -o  o  o  ooo212 

(20a) 


Mxj)MXo)  -a_r 

p-(r)/p0  =  -<“>o)(2 /xo){Tj^jro(xfy}e 


(20b) 


with 


X[  =  rt  and  x'q  =  rQ,  (i  =  1,2). 


For  small  values  of  x^>  Eq.  20b  becomes, 
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-u  r 


(22a) 


xv 

p^(r)/po  =  -(oj_/ojo) ,  x±  «  1,  (i  =  1,2) 

xo 

-™_r 

or  p^(r)  =  (K/2)polpo2(k_/.tT)c-  (22b) 

where  K  =  g/p  c/\  (23) 

o  o 


As  shown  in  Appendix  A,  Eq.  22b  is  the  asymptotic  far-field  form  of 

the  unsaturated  (pre  shock)  plane  wave  solution  of  Burgers'  equation 

27 

(Eq.  2)  derived  by  Naugol'nykh,  Soluyan  and  Khokhlov. 

Alternatively,  for  spherically  spreading  waves  centered  at 
range  r^,  since  P_(r)  =  rp^(r)/rQPo,  the  asymptotic  far-field  solution 
which  corresponds  to  Eq.  20b,  can  be  expressed  as. 


I  (X")l  (x")  -a  r 

pl<r)/po  -  -  (al_/u)o)(2/x^HT^iTTL|T)(ro/r)e  ‘  .  (24) 

For  small  values  of  x”,  Eq.  24  becomes, 

xV*  -“Ct  r 

D^(r)/po  =  -  -}(ro/r)e  “  ,  XJ  «  1,  (i  -  1,2).  (25) 


Equating  this  expression  to  the  asymptotic  far-field  form  of  the  unsaturated 

(pre  shock)  spherical  wave  solution  of  Burgers'  equation  (Eq.  1),  as 

specified  by  Eq.  A-19  of  Appendix  A,  x'J  arid  x"  become, 

i  o 


Xi  =  °iEl^aTro^’  xo  =  °oEl^aTro^’  ^  = 


(26) 


Up  to  this  point  the  discussion  has  been  confined  to  asymptotic 
far-field  solutions  of  Burgers'  equation  for  infinitely  plane  and  spherical 
wave  fields.  However,  these  particular  solutions  can  be  combined  to  obtain 
an  expression  for  the  asymptotic  far-fie!Ld  form  of  the  difference-frequency 
signal  generated  in  the  medium  by  a  dual  frequency  parametric  source  of 
finite  dimensions. 
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For  such  sources  unsaturated  (pro  shock)  solutions  of  the 

28  29 

three  dimensional  nonlinear  wave  equation  obtained  by  Muir  and  Blue 
3 

and  the  author,  show  that  if  the  parametric  interaction  takes  place 

primarily  within  the  collimation  distance  r  ,  the  difference-frequency 

2 

pressure  is  proportional  to  gj  in  contrast  to  the  dependence  on  w 
predicted  by  Eqs.  20a  and  20b.  A  heurestic  explanation  of  this  dis¬ 
crepancy,  due  to  Mellen  and  Moffett,  is  that  in  the  case  of  a  bounded 
source,  the  difference-frequencv  field,  like  the  primarv  waves,  has  a 

collimation  distance  r  which  is  shorter  than  r  on  account  of  its 

o_  o 

longer  wavelength,  i.e.,  r  =  S  /X  =  r  (u)  /w  ) .  Thus,  if  the  near- 

o  o—o  —  o 

field  primary  wave  absorption  loss  (a,^)  is  such  that  nonlinear  inter¬ 
action  occurs  primarily  within  r^,  the  difference-frequency  pressure 
field  described  by  Eq.  20b  (which  applies  to  infinite  plane  waves)  must 
be  modified  to  include  the  influence  of  spherical  soreading  losses,  so 
that  it  becomes, 


,  MxUMxU 

lp-(r)/p0l  "  fa_/M0)(2/x;Hy'  (  .)T  (  .)Hr  /r)< 

o  A1  o  K  'l 


-a  r 


i^xpi^x’) 


-a  r 


=  <«-_/%)  (2/x;)(T-Tq-)v{x,)Hro/r)e  "  ,  Vq  >  1 


(27) 


where  and  x  are  defined  by  Eq.  21. 

Alternatively,  if  the  near-field  primary  wave  absorption  loss 
(°lTro^  *s  verX  small,  most  of  the  nonlinear  interaction  takes  place 
beyond  r  ,  so  that  if  the  contribution  to  the  difference-frequency  field 
generated  within  rQ  is  neglected,  the  asymptotic  far-field  solution  is 
given  by  Eq.  20b,  which  is  reexpressed  here  for  the  benefit  of  the  reader  as. 
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P'  (r)/p 


1 , (xV) 1 .  (x")  -U  r 

c.  -  <.>0)«/,«).Ti5in^).,,r(>/r).  -  ,  Vo  <  1. 


(28) 


where  x"  and  are  defined  by  Eq .  26. 

Tn  the  general  case,  when  hotii  the  near  and  far-fields  of 

the  primary  waves  contribute  significantly  to  the  nonlinear  interaction 

process,  a  combined  form  of  Eqs.  27  and  28  can  be  expressed  as, 

i  i  n  (xQ)  -a  r 

,P_(r)/pJ  =  (-Wu)o)  (2/xo){“i— (Xi>  j  (x2>}(ro/r)e  ,  1  <  n  <  2,  (29) 

where  xi»  XQ  and  n  are  obtained  by  matciiing  the  unsaturated  forms  of  the 
asymptotic  solutions  (for  greater  than  and  less  than  unity)  defined 

bv  Eqs.  27  and  28,  as  shown  in  Appendix  B.  Thus, 


Xi 


=  o  . 


iA’X0  =  V>  d-1.2) 


i,.,r 


where 


/  = 


Fl(Vo)G  "  1/aTro’  for  Vo  >  !* 


Likewise,  n 

where  A ' 


lo8in(A/A') 

1  +  {- - 5 


log^g (^Q/w_) ' 
aTro 

El^aTro^e  1/aT^rn’  IOr  aTr'  >  !»  r' 

-*•  1  °  io  To  o 


(30) 

(31) 

(32) 


ro(wo/c *)_).  (33) 


It  is  clear  from  inspection  of  Eq.  32,  with  the  aid  of  Eqs.  31 
and  33,  that  the  index  n  approaches  unity  for  <  1,  and  approaches 

the  value  two,  for  >  1,  as  required.  The  index  n,  which  defines 

the  frequency  response  of  a  parametric  array  is  shown  in  Fig.  3  as  a 
function  of  ot^  for  different  values  of  the  "freauency  downshift  ratio" 
“«>/“- *  Likewise>  the  parameter  A,  defined  by  Eq.  31  can  easily  be  obtained 
from  Fig.  2,  which  gives  A',  as  a  function  of  oiTr\  For  small  values  of 
xi »  E<5-  29  can  be  reexpressed  with  the  aid  of  Eqs.  30  to  33,  21  and  26  as 
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P„Jl  J.  1 1  ..  gp  emu'  JI;  A  ..W 1 ! 47  ■>  4-  WWHUfly  V'M? 


-a  r 


|  p^(ir) /pQ  |  =  (“>0)n{  (a1>2/2)li1(aTro)exp(aTro)}(ro/r)c  ~  ’  x±  <K  1  ’  (i  = 

(34) 
(34a) 


-a  r 


(W-/(0r)  (°l,2/2aTro)(ro/r)e  *  aTro  1 


”  (<l)  /<*')(  .,/2)F.  (n  r  )(r  /r)e 

-  o  i  ,  2  I  I  o  a 


-a  r 


.  a  r  ] 
I  o 


(34b) 


where  Eqs.  34a  and  34b  are  the  unsaturated  forms  of  Eqs.  27  and  28, 
respectively,  as  required,  Eq .  34a  is  Vestervelt 1 s2  unsaturated  (pre 
shock)  a.svmptotic  far-field  solution  for  an  "absorption  limited"  para¬ 
metric  array  (where  the  nonlinear  interaction  takes  place  primarily 
within  the  near  field  of  the  nrimarv  waves),  and  Eq .  34b  is  the  equi¬ 
valent  pre  shock  asymptotic  solution  derived  bv  Carv8  and  the  author8,9 
for  a  "spreading  loss  limited"  parametric  array  (where  the  nonlinear 
interaction  takes  place  primarily  in  the  far-field  of  the  primary  waves) 
Equation  34  which  combines  these  unsaturated  (pre  shock)  asymptotic 
solutions  is  equivalent  to  Eq.  B-4  of  Appendix  B  as  previously  exhibited 
m  Table  Ja  and  in  Fig.  2,  where  it  was  shown  to  be  in  good  agreement 
with  the  alternative  approximations  derived  by  Mellen  and  Moffett,1,2 
Berktay  and  Leahy. ^ 

Ihe  subject  of  this  paper  however,  is  the  more  general  Eq.  29 
which  includes  the  unsaturated  (pre  shock)  solution  represented  by 
Eq.  34  as  a  special  case.  If  Xl  =  x2  =  X,  Eq .  29  becomes, 

2 


\p'Jr)\  =  (u  /cu  )n(2/Kk  /  ){_llii}2  £ 
O  o  IQ(X)  r 


ft  r 


J  K  =  8/p  c‘ 
o  o 


(35a) 


-a_r  -a  r 

(m_/wo)n(l/2Kko/)x2  1  =  (u>o)n(Kkp2)(p1>2ro)2  |  ~  ,  x 


<< 


1,2) 


(35b) 
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where  x  -  0Vli2ro)/  ;  Pj  2  -  Po|  *  |.„2. 


In  this  form,  the  amplitude  dependence  of  the  asymptotic  far- 

field  solution  is  controlled  by  the  function  {  I ^ (y) /I  (x ) }  ,  which  as 

shown  in  Fig.  4  resembles  the  well  known  characteristic  of  a  simple  two 

plate  vacuum  tube,  passing  from  a  square  law  dependence  on  the  mean 

primary  wave  amplitude  p1  9  through  a  lower  power  dependence,  to 

complete  independence  or  saturation  at  high  values  of  x(or  p.  „) . 

1 ,  z 

In  addition  to  defining  the  far-field  amplitude  and  frequency 
response  of  a  parametric  array,  Eqs.  29  and  32  can  also  be  used  to 

determine  the  far-field  directivity  function  for  the  case  of  a  "spreading 

•  m  30 

loss  limited"  arrav,  by  means  of  Lockwood's  approximation,  as 

1 3 

previouslv  utilized  hv  the  author.  This  consists  in  replacing  xi 

in  Eq.  29  by  x^KM)*  where  Di  (0 , 4> )  (i  =  1,  2) 

are  the  normalized  far-field  directivitv  functions  of  the 

individual  primary  waves,  respectively.  If  for  example, 

-3 

Xi  =  X2>  ancl  axr0  =  Np’  t^ie  dif ference-frequencv  directivity 

function  generated  in  the  medium  by  a  circular  Diston  source  with  k  a  =  10. 

o 

and  fo/f-  =  ^  is  dePicted  in  Fig.  5  as  x  varies  by  orders  of  magnitude 

3 

from  0.1  to  10  ,  respectively.  These  results  show  that  as  x  increases 


the  difference-frequency  beam  pattern  becomes  blunted  around  the  axis,  and 
the  sidelobe  levels  increase  relative  to  the  major  lobe  until  x  =  10 


when  the  beam  becomes  essentially  omnidirectional.  This  phenomenon  which 

has  been  discussed  and  confirmed  experimentally  by  Lockwood,  Muir  and 
31 

Blackstock  for  monofrequency  sources,  is  due  to  the  dependence  of 
finite-amplitude  absorption  on  the  primary  wave  amplitudes.  Since  these 


are  greatest  on-axis  in  the  far-field  of  a  directive  source,  the  amount 
of  finite-amplitude  absorption  is  also  greatest  on-axis,  and  decreases 
off-axis  as  determined  by  the  primary  wave  directivity  functions. 

In  contrast  to  the  above  case  of  a  "spreading  loss  limited" 
parametric  array  (a,j,r  <<:  1) ,  no  satisfactory  expression  for  the 
directivity  function  of  an  absorption  controlled  array  (u^r  >>  1)  as 
a  function  of  x  lias  vet  been  obtained.  However,  it  follows  from  Fig.  A 
that  it  should  he  possible  to  use  either  Merkl  inger ' s  approximation,"*2 
as  discussed  by  Childs,33  or  liartram's  series  approximation,5  for  values 


of  x  less  than  100,  where  complete  saturation  occurs. 

Having  thus  obtained  a  general  solution  for  the  asymptotic 
far-field  pressure  of  a  parametric  array,  the  next  section  will  consider 
the  construction  of  simple  scaling  laws  from  this  solution  to  facilitate 
the  design  of  parametric  arrays. 
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2.  PARAMETRIC  SCALING  LAWS 


The  asvmptotic  far-field  nressure  of  a  parametric  array,  derived 
in  the  previous  section,  can  be  expressed  in  terms  of  rms  quantities 
(denoted  by  a  bar)  by  writing  Eqs.  29  to  31  in  the  form, 


-a  r 


where 


p!(»  )/P0l  =  d//2)  (f_/f0)nF(x0,X1»X2)  (r0/t)e  ,  (37) 

(38) 


MxjH.txo) 
:'(Xo’xrx2)  =  (2/xo){y77vTTx7)' 


and 


7;  o  =  Kk  r>  r  ;  K  =  R/o  r 

n  r~\  r-\  r-\  *  ’ 


o  o  o  o  o  o 


o  o 


(39a) 


with 


*i  ■  V  =  °1  ■  KkoPoiro  *  ^  Kko"oiro; 


(39b) 


the  index  n  being  defined  by  Eqs.  32  and  33. 


Let 


SL_  =  20  log10|p'(r)  *  r  exp(u  r) 


(40a) 


and  SLq  =  20  l°gl0|poro|  with  SL.  =  20  log10 1  Po.ro |  ;  (1=1,2),  (40b) 

P2  +  2 

where  SLq  =  SL.  +  10  log1Q(— -  2  °2);  U-1,2)  (41) 


SL_,  SLq,  and  SL_^  (i-1,2)  are  thus  the  source  levels  referred  to  1  m  of 
the  difference-frequency  signal,  the  combined  primary  wave  fields,  and 
the  individual  primary  waves,  resDectivelv. 

Expressing  Eq.  37  in  logarithmic  form  with  the  aid  of  Eqs.  40 
and  41  thus  gives  the  conversion  efficiency  of  the  parametric  source 


(SL_  -  SL^)  as, 
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-  ~  . 


.......  ...... - 


(42) 


SL_  -  SLq  =  20n  log1(J(f_/lo)  I-  20  !oKl0(l-7/2) 

where  F  is  defined  by  Eq.  38. 

Since  xo  and  xA  are  required  to  evaluate  F,  it  is  convenient 
to  express  them  also  in  logarithmic  form  so  that  Eqs.  39a  and  39b  become, 

20  log10Xo  =  SLo  +  20  log10fo  +  20  loR10  +  20  lo810(N//2) 

20  log10Xi  =  SLi  +  20  log10fo  +  20  log10  +  20  logl0N;  (i=1»2) 
where  N  =  (2/2irC/P()c^)  (103) ,  for  f  in  kHz.  (44) 

If  XQ.  <,;  1  then  Eq.  42  becomes, 

SL_  *  SL0  +  20n  log(f_/fQ)  +  20  Jog10(Xlx2/2/2xo) 

=  SL1  +  SL2  +  20  log1()f_  +  20  log10A'  +  20  log1()(N/2), 

from  Eqs.  32,  34a,  and  34b,  (45) 

where  Eqs.  45  is  the  unsaturated  solution  in  logarithmic  form  previously 

3 

obtained  by  the  author. 

Since  £,  as  defined  by  Eq.  31  is  a  function  of  a„r  only,  and 

To 

N,  as  defined  by  Eq.  44  is  a  constant  for  a  given  fluid,  it  is  evident 

from  inspection  of  Eqs.  42,  43a,  and  43b  that  by  scaling  the  source  levels 

to  1  kHz,  according  to  the  procedure  introduced  bv  Mellen,  Konrad,  and 

Browning,34  the  conversion  efficiency  (SL  -  SL  )  can  be  redefined  as  a 

o 

function  of  the  scaled  source  level  SL*  and  the  nondimensional  parameters 
aTro  and  fQ/f_.  Thus  if 

SL*  =  SL_  +  20  log1Qfo,  fQ  in  kHz  (46a) 


(43a) 

(43b) 
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and  if 


SL*  =  SL  +20  log.nt  ,  SL*  =  SL.  +  20  log._f  ,  (1=1,2)  (46b) 

o  o  10  oi  i  ft10  o  ’ 

then  Eqs .  42,  43a,  and  43b  become, 

SL*  -  SL*  =  SL_  -  SLq  =  2 On  log^^/O  +  20  log1(J(F//2)  (47) 

20  log1QXo  =  SL*  +  20  log1()7  +  20  1or1()(N//2)  (48a) 

20  log1QXi  =  SL*  +  20  log1()A+  20  LoBj(JN,  (1=1,2).  (48b) 

Curves  of  the  conversion  efficiency  (SL  --  SL  )  and  SL*  obtained  from 

-  o  - 

Eqs.  38,  47  and  48  with  SL.  equal  SL„ ,  are  shown  as  functions  of  SL*  for 

i  z  o 

particular  values  of  ffl/f_  and  a.r  in  Figs.  6a,  6b,  7a,  7b,  8a,  and  8b, 

where  the  medium  is  water  (20  log^  N  =  281  dB//luPa)  ,  all  pressures 

being  referred  to  1  pPa.  Using  these  scaled  characteristics,  parametric 

sonars  with  similar  values  of  fQ/f_  and  a  r  ,  but  with  different  operating 

frequencies  and  source  dimensions  can  be  compared  as  functions  of  the 

scaled  source  level  SL*. 

o 

At  this  point  it  is  instructive  to  compare  the  conversion 

efficiency  defined  by  Eqs.  38,  47,  and  48  with  results  obtained  from 

1  2 

Mellen  and  Moffett's  model.  ’  Such  a  comparison  is  shown  in  Fig.  9 

where  fQ/f_  =  10  and  varies  from  1  to  10  Np.  It  should  be  noted 

that,  unlike  Figs.  6b,  7b,  and  8b,  the  ordinate  in  Fig.  9  gives  the 

conversion  efficiency  as  (SL_  -  SL^  ^ )  which  is  referred  to  the  mean 

primary  wave  source  level  SL1  2>  in  keeping  with  Mellen  and  Moffett's1’2 

presentation  of  their  results.  Likewise,  the  abscissa  in  Fig.  9  is 

given  in  terms  of  the  scaled  mean  primary  wave  level  SL*  It  follows 

I ,  z 

therefore,  from  the  discussion  at  the  beginning  of  this  paper  that  the 
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Fig.  7a 
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SLq,  Scaled  Primary  Source  Level 
(dB  //  1m Pam  kHz;  rms) 


difference  between  the  results  of  the  two  models  shown  In  Fig.  9,  which 
becomes  more  accute  as  SL*  ^  increases,  can  he  attributed  entirely  to 
the  higher  order  spectral  interactions  neglected  in  Mellen  and  Moffett's 

1  2 

model.  '  These  interactions,  which  are  taken  into  account  in  the  model 

defined  by  Eqs.  38,  47,  and  48,  give  rise  to  enhanced  finite-amplitude 

losses,  thus  causing  the  conversion  efficiency  to  decrease  more  rapidly 

(after  its  maximum  value  has  been  readied)  as  the  primarv  wave  source 

levels  are  increased.  Numerical  solutions  of  the  spherical  form  of 

Burgers'  equation  (Eq.  2,  £  =  1)  obtained  with  the  aid  of  a  comouter 

9 

program  previously  developed  by  the  author  are  also  shown  in  Fig.  9 
_2 

for  aTrQ  <:  10  Np,  which  defines  the  range  of  validity  of  the  spherical 
wave  model.  Although  extremely  unstable  and  difficult  to  generate  at 
high  scaled  source  levels,  these  numerical  solutions  tend  to  confirm 
ttie  trend  of  the  analytical  far-field  model. 

In  order  to  further  simplify  the  presentation  of  scaled  results, 
a  higher  level  of  scaling  can  be  introduced.  Thus, 


and  if 


SL**  -  SL*  +  20  log1Q/  -  2 On  log1()(f_/fo) 

=  SL_  +  20  log1()fo  +  20  log1()7  -  20n  log10(f_/fQ)  (49a) 
SL**  =  SL*  +  20  log10^ 


then  Eqs.  47  and  48  become, 

SL**  -  SL**  =  (SL_  -  SLq)  +  20n  log^f^ij  =  20  log10(F//2)  (50) 

20  1o£:10*o  =  SLo*  +  20  ^10(N//2)  (51a) 

20  loR10Xi  =  SLf  +  20  lo^i0N»  (i-1* 2) .  (51b) 


34 


Ising  Eqs.  38,  50,  5La,  and  51b,  (SI,**  -  Si,**)  and  SI,**  can  be  evaluated 
as  functions  of  SI,**  to  give  the  "universal"  scaled  characteristics  shown 
in  Figs.  10a  and  10b.  Experimental  data  obtained  by  a  number  of  different 
investigators34  38  under  entirely  different  conditions  with  a  variety 
of  parametric  sources,  has  been  scaled  and  superimposed  on  these  character¬ 
istics,  showing  reasonably  good  agreement  despite  the  scatter.  A  summary 
of  the  source  parameters  used  to  obtain  this  data  is  given  in  Table  III, 
which  also  includes  values  of  20  log1Q(/)  calculated  for  the  benefit  of 
the  reader  from  F.q.  33.  With  the  exception  of  experiments  1,  2,  3,  and 

7  of  Table  III,  the  nearfield  primary  wave  absorption  loss  2a  r  was 

oo 

given  or  could  be  calculated  from  reported  values  of  a  and  r  .  Since 

o  o 


experiments  1,  2,  and  3  were  performed  in  a  hot  salt  water  medium, 
Schulkin  and  Marsh' s3'  expression  for  acoustic  attenuation  in  sea  water 
was  used  to  calculate  the  values  of  aQ  shorn  in  Table  III,  based  on  an 
assumed  temperature  of  80°F  with  a  salinity  of  30  parts  per  thousand. 

In  order  to  check  the  accuracy  of  these  assumed  conditions,  the  attenua¬ 
tion  coefficient  for  experiment  #4  (which  was  also  conducted  in  the  same 
salt  water  environment),  was  calculated  using  the  same  values  of  tempera 
ture  and  salinitv  to  give  the  value  1 . 17  (10_2)N’p/m  shown  in  Table  III. 
This  calculation  is  in  very  good  agreement  with  the  value  of  1.2(10_2) 
Np/m  quoted  by  the  investigators,35  who  also  gave  the  attenuation 
coefficient  for  the  250  kHz  source  in  experiment  //2  as  10_2Np/m  in 
keeping  with  the  calculated  value  of  0. 95 (10‘2)Np/m  shown  in  Table  III. 

In  the  case  of  experiment  #7,  since  the  medium  was  fresh  water  at  normal 
temperature  and  pressure,  a/f 2  was  assumed  to  be  25(l0-15)Np  secsV 
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TABU'  I  I  I 


t 

w* 

9 

Experiment 

legend 

f 

o 

kHz 

a 

o 

Np/m 

r 

o 

m 

2a  r 
o  o 

Np 

20  log1Q 

dB 

** 

1 

34 

NUSC 

• 

175 

6.76(10~3) 

16 

2 . 16 (10_1) 

3.0 

*"■ 

2 

34 

NUSC 

A 

250 

0. 95 (10-2) 

8 

1 . 52 (10_1) 

4.5 

««> 

3 

34 

NUSC 

D 

720 

2.08(10~2) 

4 

1.66 (10_1) 

4.0 

r 

i  * 

4 

35 

NUSC 

0 

330 

1. J7(10-2) 

0.25 

6.0(10-3) 

13.0 

*  ■» 

3 

NUSC36 

O 

65 

2 . 23 (] 0~3) 

28 

1 . 3 (10_1) 

5.0 

6 

37 

Muir  and  Gillette 

0 

450 

6. 25 (10~3 ) 

1.4 

1. 7  5 ( 10_ 2 ) 

11.0 

-* 

7 

A.  F.ller38 

J 

1435 

5. 13(10  2) 

0.3 

3 . 09 (10~3) 

14.0 
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An  example  of  the  calculations  required  to  scale  the  data 

shown  in  Figs.  10a  and  10b  is  given  in  Tables  IVa  and  IVb  for  the  case  of 

the  250  kHz  source  listed  as  experiment^  in  Table  Ill.  The  values  of 

34 

f  given  in  Table  IVa  correspond  to  those  of  the  experiment,  the 

values  of  20  log/'1  being  calculated  from  Kq.  33  with  a  rQ  (-  2aQr  )  as 

giver  in  Table  111.  The  index  n  is  then  calculated  from  Eq.  32.  In 

Table  IVb  the  actual  measured  values  of  SL  and  SI.  (for  two  different 

o 

difference-frequencies)  are  reproduced,  where  the  pressures  are  referenced 

34 

to  1  yPa  rather  than  1  ybar,  as  used  bv  the  investigators.  These 

measurements  are  then  scaled  using  Eqs.  46a,  46b,  49a,  49b,  and  Table  IVa. 

Finally,  the  scaled  difference-frequency  levels  are  averaged  to  give 

<SL**>,  on  the  grounds  that  the  differences  between  them  are  caused  by 

experimental  inconsistencies  which  should  disappear  on  the  average. 

Example:  In  order  to  use  Figs.  10a  and  10b  to  evaluate  the  effectiveness 

of  a  parametric  arrav,  consider  the  case  of  the  65  kHz  source  listed 

as  experiment  5  in  Table  III,  with  2a  r  =  a„r  =  0.13  Np,  and  20  log,- A  = 

oo  1  o  10 

5  dB. 

36 

If,  as  was  the  case  in  the  experiment  SL^  =  SL^  =  249  dB//l  yPam 
rms,  then  SLq  =  251  dB//l  yPam;  rms,  and  from  Eqs.  46b  and  4>o, 

SL**  =  SL  +20  log  f  -  20  log.n  ;  f  in  kHz 
o  o  10  o  iu  o 

=251+36+5 

=  290  dB. 
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TABLE  lVa 
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I 

f 

£  /f 
o 

20  1°8io(/' 

’)  n 

20n  log10(fo/f_)  20 

los10( (£o/£ 

)n(f  t)} 
0 

T 

kHz 

~ 

dB 

- 

dB 

dB 

J- 

12 

20.8 

-J  2 

1.63 

43 

95.5 

*<0» 

i 

a* 

6 

41.7 

-17 

1.66 

54 

106.5 

TABLE  IVb 
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Experimental 
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Results 

Scaled 

Exnerimental  Results 

SL  SL 
o 

(12  kHz) 

SL_ (6  kHz) 

SL**  SI, 
o 

(12  kHz) 

SL** (6  kHz) 

<SL**>  <SL**>-SL** 

-  -  o 

1. 

dB 

dB 

dB 

dB 

dB 

dB 

dB 

dB 

240 

189 

178 

292.5 

284.5 

284.5 

284.5 

-8.0 

<•  w 
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286.5 

276.4 

2/6.4 

276.4 

-10.1 

i 

228 

172 

160 

280.5 

267.4 
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266.9 

-13.6 

♦  nr 

222 

160 

150 

274.5 

255.4 

256.4 

255.9 
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Thus,  from  Fig.  10a,  SL**  =  280.5  dB.  If,  as  in  the  experiment  ,  f_  = 

3.5  kHz,  then  from  Eq.  33,  20  log/'’  =  -9  dB;  hence  from  Eq .  32,  n  =  1.55, 
with  20n  log^Q(fQ/f_)  =  39.5  dB.  Finally,  from  Kqs.  46a  and  49a, 

SL_  =  SL**  -  (20  log1Qfo  +  20  log1()A  +  20  log1Q(f  Jt _) )  ,  or 

SL_(3. 5  kHz)  =  280.5  -  (36  +  5  +  39.5) 

=  200  dB,  which  is  equal  to  the  level  obtained 
36 

experimentally.  A  very  interesting  result  which  follows  directly  from 
Eq.  50  and  Fig.  10b  is  that  the  maximum  conversion  efficiency  of  a 
parametric  source  for  a  given  frequency  downshift  ratio  (fQ/f_),  can  be 
expressed  as, 

(SL  -  SL  )  -  -  9  -  20n  log  At Ji  )  (52) 

-  o  max  10  o  - 

where  n,  as  defined  by  Eqs.  32  and  33,  is  a  function  of  a„r  and  a„r'. 

1  o  1  o 

For  the  case  of  the  65  kHz  source  previously  considered,  the  difference- 

frequency  is  3.5  kHz  and  thus,  20n  log^Q(fo/f_)  =  39.5  dB,  so  that  the 

maximum  conversion  efficiency  in  this  instance  is  -  48.5  dB.  Since  the 

36 

conversion  efficiency  obtained  experimentally  was  -51  dB  for  a  source 
level  of  251  dB,  it  follows  that  this  sonar  was  operating  close  to  its 
limit  of  maximum  conversion  efficiency. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


A  new  analytical  model  has  been  derived  in  this  paper  for 
calculating  the  asymptotic  far-field  pressure  and  conversion  efficiency 
of  a  parametric  array.  Unlike  previous  models^  ^  this  new  model  includes 
the  effect  of  all  degenerative  spectral  interactions  on  the  conversion 
efficiencv  of  the  parametric  mixing  process.  Tv;o  levels  of  scaled 
performance  characteristics  have  been  derived  from  this  analvsis,  which 
should  simplify  the  problem  of  designing  small  scale  model  tank  experi¬ 
ments  to  simulate  the  performance  (and  thus  assist  the  design)  of  large 
scale  Darametric  sonars. 


The  model  has  been  shown  to  give  good  agreement  with  experi¬ 


mental  data  obtained  by  different  investigators^  over  a  wide  range 


of  experimental  conditions.  However,  a  word  of  caution  must  be  added  at 
this  point  concerning  the  comparison  of  data  obtained  in  a  dispersive 
medium  such  as  sea  water,  with  the  predictions  of  analytical  models 
derived  from  a  nondispersive  wave  equation. 


Unfortunately,  most  of  the  high  amplitude  data  available  in 
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the  literature  ’  which  was  used  to  test  the  model  derived  in  this 


1-5 


paper,  and  all  previous  models,  was  obtained  in  a  salt  water  medi  m 
at  primary  wave  frequencies  between  175  kHz  and  720  kHz.  Since  this 
range  of  frequencies  lies  within  the  'resonant'  region  of  the  salt 
water  relaxation  process,  the  attenuation  coefficient  is  no  longer 
proportional  to  the  frequency  squared  throughout  the  spectrum,  as 
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assumed  by  the  nondispersive  models.  Thus,  the  spectral  components 
generated  In  sea  water  by  a  parametric  sonar  operating  within  this 
'resonant'  relaxation  frequency  range  will  be  subject  to  entirely 
different  rates  of  absorption  depending  on  where  the  mean  primary  wave 
frequency  is  centered;  a  process  which  may  enhance  or  degrade  the  con¬ 
version  efficiency  relative  to  its  potential  value  in  fresh  water. 

Obviously  these  remarks  do  not  apply  significantly  to 
unsaturated  parametric  sonars  operating  at  pre-primary  wave  shock  levels 
in  sea  water,  nor  do  they  apply  to  high  power  parametric  sonars  in  sea 
water  whose  primary  wave  frequencies  are  such  that  their  fourth  or  fifth 
self-generated  harmonics  do  not  excede  100  kHz.  Likewise,  high  power 
parametric  sonars  operating  above  1  MHz  in  sea  water  may  be  unaffected 
unless  the  difference-frequency  component  and  its  harmonics  lie  in  the 
'resonant'  relaxation  frequency  range  between  100  kHz  and  1  MHz. 

However,  in  order  to  provide  effective  tests  of  the  saturation 
limits  predicted  in  this  paper,  and  compare  them  with  previously  derived 
predictions,1  5  saturated  parametric  experiments  in  nondispersive  fluids, 
such  as  fresh  water,  are  required.  Likewise,  in  order  to  predict  the 
saturation  limit  of  parametric  sonars  operating  in  sea  water  whose  spectrum 
falls  significantly  within  the  'resonant'  relaxation  frequency  range,  it 
will  be  necessary  to  include  the  effect  of  'relaxation  absorption' , 
wherever  possible,  in  the  nondispersive  analytical  models. 
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APPENDIX  A 


If  the  primary  wave  amplitudes  of  a  parametric  source  are 
significantly  lower  than  their  respective  shock  thresholds,  the  time 
waveform  will  be  relatively  undistorted.  Under  these  conditions  it  is 
possible  to  solve  Eq.  2  by  means  of  the  method  of  successive  approxima¬ 
tion.  Basically,  this  method  consists  in  solving  the  linear  part  of 
the  equation  by  neglecting  the  nonlinear  term  to  obtain  the  "first 
approximation,"  which  is  the  undistorted  time  waveform.  Substituting 
this  solution  in  the  nonlinear  term  and  integrating  the  resulting  inhomo¬ 
geneous  equation  thus  gives  the  "second  approximation."  The  process  can 
then  be  repeated  to  obtain  higher  approximations. 

In  or  ar  to  simplify  the  procedure  of  selecting  the  nonlinearly 
generated  spectral  components  of  interest  it  is  convenient  to  take  the 
Fourier  transform  of  Eq.  2  so  that  it  becomes. 


where 


and 


df+  (“2/A)l  'I3tA{p2) 

'  CO 

yP(C.t’)}  =  i*u(c)  =  P(c,t’)e"ja)t'  dt' 

—  00 

Cl'W’  ■  Xc.f)  =i-  ['i>  (Oe^1'  d„. 

—  00 


(A-l) 

(A-2) 

(A-3) 


Solving  the  linear  part  of  Eq.  A-l  for  a  dual  frequency  parametric  source 
operating  simultaneously  at  angular  frequencies  u>1  and  ^  with  peak 
amplitudes  Pq1  and  Pq2  respectively,  and  using  Eq.  A-3,  gives  the 
undistorted  time  waveform. 


It  should  be  noted  that  Eq.  A-10  is  the  approximate  solution  originally 
obtained  by  NaugDl'nykh,  Soltiyan  and  Khokhlov.^ 

Likewise,  for  spherical  waves  (»L  =  1)  Table  II  gives, 

4  =  oUn(r/ro),  so  that  d;  =  o^(dr/r)  (A-13) 
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APPENDIX  B 


In  order  to  determine  the  parameters  X^»  X2  an<*  n  so  that  Eq.  29 
will  approach  Eqs.  27  and  28  for  the  asymptotic  extremes  of  aTrQ  greater 
and  less  than  unity,  respectively,  it  is  expedient  to  begin  with  the  form 
of  Eq.  29  for  small  values  of  ^  which  is. 


,  ,  „  X1X1 

I  Pl<r)P0  I  =  (<*>_/%,)  — Hro/r)e  »Xi<<  1,  (i  ■  1,2).  (B-l) 

o 

Likewise,  for  small  values  of  x|  (as  defined  by  Eq.  21), 

Eq.  27  can  be  expressed  as, 


p^(xr ) /pQ|  =  (w_/ojo)  (  (a1  2/2) 


r  -a  r'  ,  -a  r 

e  — }(ro/r)e  ,  aTr  »  1  (B-2) 

o  o 


->  (w  /a)  )  ( o  ,/2a  r  )(r  /r)e 
— o  1,2  to  o 


-a  r 


Again,  for  small  values  of  xV  (as  defined  by  Eq.  25),  Eq.  28 
can  be  expressed  as, 


p '  (r ) / p  |  =  (w  /oj  ){(o.  9/2) 
—  o1  —  o  1,4 


"aTr'  dr'  ^-r 

e  —H  (r  /r)e  ,  aTr  «  1  (B-3) 

r  r  o  i  o 

o 

-a  r 


(a>  /oj  )(a.  „/2)E  (aTr  )(r  /r)e 
-o  1,4  ilo  o 


Combining  Eqs.  B2  and  B3  thus  gives, 


p^(r)/pj  =  ((0_/a)o)(a1  2 / 2 ) {  (aj  L  ) 

’  "  °  '  o, 

f  OD  ^Tr 


rro  ^Tr’  dr'  ^  r  "V  dr- 
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77— r  rr  dr'}(r  /r)e  ;  r'  =  r  (00  /a)  ) 
rrr  O  000- 
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Equation  B-4  which  was  previously  derived  by  the  author  as  shown  in 
Table  la  can  be  reexpressed  in  terms  of  as  follows: 


|p'(r)/p  |  =  (ui  /  to  )n(o  _){E  (u  r  )exp(a  r  ) }  (r  /r)e 

~  O  —  O  1,4  JL  L  O  iO  O 


(B-5) 


where  the  index  n  is  obtained  by  equating  Eqs.  B-4  and  B-5  to  give. 


JLog  (A/A') 

n  =  1  +  — -} 

ioglO^o/u).) 


(B-6) 


with  A  =  E  (a  r  )exp(a  r  )  and  A'  =  E  (a  r ' ) exp(a  r ' ) , 

JLIO  iO  ilO  lO 


(B-7) 


Equating  Eqs.  B-l  and  B-5  thus  gives  and  XQ  as. 


Xi  =  °iA  and  -  V 


(B-8) 
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